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ABSTRACT 



A ship in steering is considered as a physical mechanism that is 
forced by a rudder movement to produce a response. Emphasis is 
laid upon the relation between the forcing and the response, leaving 
aside any detailed consideration of the forces concerned. 

Course-Keeping with Automatic Control techniques is studied 
following this concept. Computer programs are developed to simulate 
different conditions. Interpretation of the results is made to evaluate 
the different methods used. 
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I. INTRODUCTION 



Problems that justify the use of Automatic Control for ships: 

The maneuverability of a ship is determined by the physical 
properties of the ship in its environment and the physiological and 
psychological properties of its navigator. 

A successful mission of a vehicle in its operating environment is 
a function of its handling as well as its design. The relationship 
between helmsmanship and the steering system has been discussed in 
detail by Abkowitz [1], Brard [2], Davidson [3], Nomoto [4], among 
others, and the conclusions indicate that the hydrodynamic design 
factors alone can not significantly improve the handling qualities of 
ships . 

The past decade will be known as the decade of the Super- Tanker s . 
Ships were made larger for many reasons, one of which was to im- 
prove operations efficiency. To realize the full benefit of the larger 
size, the ship had to be steered as efficiently as possible. Almost 
every Super- Tanker built nowadays is course unstable and the control 
of such a system has been of great concern. We can see easily that 
it is not difficult to design an autopilot if enough rate control is 
available [ 5 ] . 

The basis for the need of Automatic Control was a minimum 
number of operator controls, minimum rudder orders for course- 
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keeping and automatic course correction without overslioot. But the 
main reason for this type of control is that steering an unstable ship 
by manual control requires a helmsman to detect a very small amount 
of rate of turn and to react to this in time with the proper correction. 

Even though the study of man as a pilot of giant tankers is rel- 
atively young, yet it is evident that any pilot, not withstanding his 
experience and seamanship, has some physiological and psychological 
limitations relevant for the execution of his task. 

Several studies carried out [6 and 7] provided us with some in- 
sight. First a purely physiological problem is a limited capacity for 
the perception of motion. Very slow motions cannot be perceived by 
human beings as is the case with the hands of a clock. To be more 
specific, the smallest yaw velocity human beings can perceive is 
about 1 minute of arc per second. The smallest just noticeable 
acceleration or deceleration of the yaw motion occurs when the yaw 
velocity is doubled or halved within 5 seconds. Analysis of ship man- 
euverability experiments on full-scale models revealed that tankers 
above 100, 000 Ton. move so slowly during maneuvers that many of 
the accompanying velocities and accelerations are not perceptible 
for the man on the bridge. This means that some information which 
is indispensable for the correct execution of the maneuvers is not 
available for the navigator due to his own physiological limitations. 
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The inertia o f large tankers also raises a problem of a more 



psychological nature. It is well known to pilots and Captains that 
there is a time lag between a rudder command and the reaction of 
the ship. This lag can be easily of the order of 15 seconds. Psy- 
chological research on tracking behaviour revealed that steering 
tasks raise difficulties when the time lag is longer than 4 seconds. 
Also related to inertia is the ability to anticipate motions over long 
time intervals. Even very simple maneuvers take a long time and 
should therefore be initiated a long time beforehand. It is however, 
very difficult to anticipate maneuvers over long time intervals as 
human beings have only a limited capability for extrapolation of slow 
motions, i. e. , the navigator can hardly predict what his position 
will be in 10 minutes and therefore is not able to initiate the neces- 
sary correction at the required moment. 

Still another factor is the ability to discriminate among the 
different components of the ship's motion, those due to wind, current 
and steering movements initiated by the pilot. 

From experiences in working with the concept of the time con- 
stant of a ship, conclusions resulted in that the difficulty in steering 
depends largely on this parameter, since for the same statistical 
characteristics, the unstable ship with the small value of absolute 
time constant, changes its course loo quickly to be followed by the 
human operator. 
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II. COURSE- KEEPING 



A. DEFINITION 

Course-Keeping quality is defined as the ability of a ship to keep 
course easily on a seaway. This is not identical to stability on course, 
because the former depends upon the behaviour of the steering devices 
as well as of the ship, while the latter is a thoroughly passive char- 
acter of the ship. So the problem of course-keeping should be treated 
by analyzing a closed-loop feedback system composed of a ship, an 
automatic or manual steering control unit and steering gear. It is 
possible, however, to show simply that good stability on course 
yields a good course-keeping quality, in general, in this form: 

1. A ship that is more stable on course has less occasion to 
deviate significantly from its course because the stimulated motion 
decays more quickly even with the rudder amidship. 

2. A ship that is more sta.ble on course responds more quickly 
to steering, so that any coxirse deviation can be corrected more 
easily. 

The stability on course is related to the decay of the yaw and 
sway of the ship with a rudder amidship after being disturbed for a 
short while. In the case of an unstable ship this results in a negative 
damping of yaw and sway. Stability on course is also related to quick 
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response to steering. This is of particular importance in considering 
the cours e- keeping quality of a ship. 

B. STEERING QUALITIES OF SHIPS 

The dynamic stability on course of an unsteered ship is measured 
by the number of ship lengths traveled by a stable ship, in the time 
required to reduce an accidental deviation from undisturbed motion 
in a straight line to 1/e of its initial value, the rudder remaining in 
the a mid ship position, this quantity is defined as l/|p^J . 

A negative value of pj means that the ship is dynamically stable, 
i. e. , when a dynamically stable ship moving on straight course is 
disturbed slightly, it settles down on a new straight course that is 
close to the original one (Figure 1). The greater the negative magni- 
tude of p, , the more rapidly the ship settles on its new course and 
the closer the new course is to the original one. A positive value of 
p^ means that the ship is dynamically unstable, i. e. , when a dynam- 
ically unstable ship moving on straight course is disturbed slightly, 
it will go into a steady circuling motion, even though the rudder is 
held amidship. 

The primary piirpose of an automatic steering device, just as of 
a human helmsman, is to compensate for disturbances, making the 
ship maintain a prescribed heading and thus making it directionally 
stable. 
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FIGURE 1 



C. UNSTABLE HULL 

In this study of Automatic Course-Keeping we will he working 
with a 200, 000 DWT super- tanker of the following characteristics: 
Length 310 meters 

Breadth 47. 16 meters 

Draft 18. 90 meters 

Steering Quality indices: 

= -269. 3 seconds 

= 9. 3 seconds 

= 20.0 seconds 

K = -0. 0434 rad/sec 

Maximum Rudder deflection 30° 

Maximum Rudder rate 2.32 degrees /s econd 
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D. AUTOMATIC CONTROL AND COURSE-KEEPING 



One of the functions of ship control is to maintain a ship's 
heading. In performing this function, a helmsman deflects the rudder 
in a way which will reduce the error between the actual and desired 
heading (Figure 2) 



A good helmsman will not only deflect the rudder in response to 
the heading error, but he is also sensitive to the angular velocity of 
the ship and he will ease off the rudder and apply a little opposite 
rudder in order to prevent overshooting the desired heading. It 
follows that an automatic control should also be responsive to control 
signals measuring both error and angular velocity. Thus a rudder 
under automatic control might be deflected in accordance with the 
following linear expression: 







Actual 

UCAOttOC 



FIGURE 2 
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where 



& ~ s -t- <f 

>Jji yy 

£ helmsman's rudder angle 

angular velocity of ship 
'yy heading error 

The basic action of the automatic control or "autopilot 11 is called 
proportional control which means to give a helm angle proportional 
to the amount of course deviation and sometimes also proportional 
to the time rate of the deviation. 

In general the complete steering system can be represented as 
in figure 3 
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FIGURE 3 



The response of a ship to steering or, the maneuverability of a 
ship is usually described by a set of equations of motion of side- 
drifting velocity and turning-angular rotation: 
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It is however, more convenient to use a single equation of motion 
describing immediately a relation between turning angular motion and 
steering, because control signals in autopiloting relates only to turn- 
ing angular motion. 
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III. TRANSFER FUNCTION IN STEERING 



The general principle of linearizing the equation of motion of a 
ship in steering may be summarized as follows: 

1. The coupling of surging to sway and yaw is ignored because 
of its second-order effect on the latter motions [ 19 ]. 

2. Hydrodynamic forces acting upon the ship are expressed as 

linear functions of angle of drift, i. e. , sway velocity, curvature of 

the ship’s path, i. e. , yaw angular velocity, and angle of helm. 

* 

3. Hydro-inertial forces are expressed as linear functions of 
the acceleration and angular acceleration of the ship. 

Linear equations of motion of this kind have been introduced in 
a number of articles [8, 9, 10 and 17], and have proved promising in 
interpreting and predicting behaviours of ships in steering. In the 
notation used here, they arc written as: 



( L A) ( yY) ' + tof) h -h Yfi P - (Yer) (yy}' 



Y'r)e*yJs 




F*) & ‘ L ( L /v ) A/(3 & - 



A/'fi v = a/j cf 



III- 1 



Using these equations we can obtain the yawing (0 (t)) and the 
swaying (/* (t)) in response to any movement of the rudder ( £ (t)). 

Nov/ limiting our interest to the steering- to-yaw response which, 
in general, is naturally of major importance, we introduce the 
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transfer function in steering. Taking the LaPlace transform of both 



sides of the equations III- 1 , we obtain: 

- ( yW ( Vy ) P ( °) + L ( L /ir) ( to 7 \My ) 5 + y? ] P 
( to -t to'* - Yy ) © ( S) — Ys c5 c 5 ) 

- ( y<rf( Z zz + J ?z) ® (o ) +[( i *r)'(T' gg +Jt r )s +(%/<]ecs) 

- ///? s) - a/j J ( sj 

• * 

where j}> {0) and Q (o) denote p and 0 respectively at t=0, 
namely at the beginning of a rudder movement, and where 

ft ( 5 *) is the LaPlace transform of /3 

0 

Q (S) is the LaPlace transform of © 

^ (S ) is the LaPlace transform of cf 
Taking into account the following relation at t = 0 

( L /f) (Z'zz ®t°) +■ ( L /tr) Aj' r 9 co) - a// 1 fi(o) = o 

we can eliminate and {o) from the transformed equations [ 19 ], 

to yield a single equation that describes the steering- to -yaw re- 
sponse, that is: 

£(S) +. s 1 a T ± 0(O ) m ?_ 

c H-T, s )( nns) ( /+r, s ; ( t + T t s) 



Since later in this study we are going to be dealing with auto- 
matic piloting, and the control signal under this condition relates 
only to turning angular motion. 
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where 



l< _ / V/ L ) N r- Yj ■+ Yfi M's 

Yp A >r - ( YP H P 

T ~ t ~r L __ ) C — ~ — 1 1 r iiSPjll. ~!d j. tsL 'l i * 

Ypf's’f ~ ( TT}' L ~Y r) ^ f* 

T, T.Jl/vY .JPl^jdPlidhl 

( Ka' -f- \<A-'y) dj 

v* yj * w 

The first term of the right side of equation III- 2 corresponds to a 
motion excited by steering and the second to another resulting from 
the initial motion at t=0, the latter disappears if the ship is running 
straight at t=0 ‘without any yaw or sway. Thus the rational function: 

&(S) _ Kd+'T+s) 

JfS) ~ ( H-T l 5)C^ T 7.s) 

describes the response behaviour of a ship in terms of the LaPlace 
transform. We may call it the transfer function of the ship in steering 
Retransforming equation III- 2 , we obtain a single differential 
equation that describes ship motion just as does the original equation 
of motion (Equations III - 1 ) as far as the steering- to-yaw response is 
concerned, that is: 

77 7i 6 (7, tTi) d_ o + o = !< S + kT s £I_ 

d t x 
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better known as Nomoto's Equation, this describes directly the yaw 
response of a ship. The four coefficients K, Tp T^, and T^, which 
are composed of the coefficients of the original equation of motion, 
constitute a set of characteristic figures representing the response 
behaviour of the ship, they are called the steering quality indices. 

A. STABILITY CRITERION OF AUTOPILOTING BY KOCHENBURGER'S 

METHOD 

A control loop describing autopiloting of a ship is illustrated in 
figure 3 in the form of a block diagram. Since the stability criterion 
of Kochenburger [11] is based upon whether sinusoidal signals grow or 
decay in circulating through the loop, it is necessary first to obtain 
the response of each element composing the system, to sinusoidal 
signals. 

The response of a ship to a sinusoidal signal (that is, in this case, 
to put a rudder sinusoidally to both sides with a certain frequeneyfw) 
may be determined through the equation of Nomoto: 



£ (t) . r A) a to L 



where 




A jj/j ) ~ C 1 * to 1 3 ; ) 



C / -f t iAj 7J ) C ( + c to 7V ) 




/? A A ( * 1 °° S 3 J 

( / + CcoT, ) c / + <~U) TO 
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and where is the amplitude of sinusoidal steering. A( uJ) is called 



an amplitude ratio and 0 (u>) a phase difference. Both of them are 

functions of frequency w only, as a common feature of linear systems 
Next, the response of an electrohydraulic steering gear, which is 
widely used for most present ships, may be described by the following 
equation: 



where T £ is a time constant of the steering gear. 

Then we get a description of the response of a steering gear to a 



Finally, we get easily the response character of a compass as 
follows, because it may be considered a simple integrating element 
transforming 0 into © : 





y t 



sinusoidal 




as follows: 




where 




A £ (^) 
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Since the transmitting character through all linear elements has 
been thus obtained, if a similar character of the remaining element, 
viz, an autopilot, is defined, we can judge whether sinusoidal signals 
grow or decay in circulating through the control loop. It is impossible, 



tial equation and then to obtain its response character in the foregoing 
manner, if considering such a discontinuous element as a weather- 



order to avoid frequent steering in practical application for rough 
seas. 

Fortunately, however, the response itself of an autopilot with a 
weather-adjust mechanism to a sinusoidal input signal (that is, in this 
case a course deviation) may be obtained easily by expanding the re- 
sponse into a Fourier series with a fundamental frequency that is the 
same of the input signal, as follows: 



however, to describe the response of an autopilot by any linear differ en- 



adjust mechanism, which is used to modified proportional control in 
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is the proportionality constant connecting a 
course deviation to a helm angle to be called 
for, 

is the amplitude of course deviation, and 

are the coefficients of Fourier series which 
may be obtained through the usual procedure 
of Fourier expansion if the form of <5 is given. 

Considering here that all linear elements are much more insensi- 
tive to a signal with the higher frequency, we can neglect all the higher 
frequency terms. This is the approximation of Kochenburger , and its 
validity depends upon how much the higher frequency signal decays 
through the linear elements and also how much higher frequency com- 
ponents are included in the original S ; viz, how much the original £ 
resembles a pure sinusoidal form. In the present case, Kochenburger ’s 
approximation may be fairly valid because a ship is quite insensitive to 
a high-frequency steering because of her large inertia, and also the 
basic action of an autopilot, even with a weather-adjust mechanism, is 
proportional control that produces a sinusoidal output in response to 
a sinusoidal input. 

B. FIRST ATTEMPT TO STABILIZE THE SYSTEM 

The autopilot, the steering engine and the ship with its rudder all 
form different components of a closed loop system, each component 
characterized by its transfer function or the complex ratio of output to 
input. The theory for such control systems and their stability has been 
developed in electric network engineering, and it is natural that the 



where C 

1 






a . a , a , and a , 

1 2’ 3 4 
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dynamic problem lends itself to studies in analog computers, where 
each component is represented by its equivalent electric circuit. The 
stability of the closed loop system may be judged from the total open 
loop response recorded at several frequencies, without a knowledge 
of the individual transfer functions. 

The conditions of directional stability in automatic steering along 
a fixed course had been discussed by Minorsky [12] by means of the 
technique of added derivatives, applied to a simplified one-degree 
of freedom oscillation and including several types of position and rate 
control. In 1946 Davidson and Schiff [ 1 3 ] took a large step towards a 
better understanding of the interrelation between the performance of 
a ship on a straight course and in turning, pointing out the nonlineari- 
ties in the behaviour of the unstable ship and including an interesting 
treatment of the transients when entering a steady turn. They also 
used linear theory for establishing a formula for the radius of the 
turning circle. 

The work of Davidson and Schiff did much to stimulate other 
authors. Among many others, a paper by Williams [14] on initial 
stage motion and a report by Schiff and Gimprich [15], who studied an 
automatic control system where the rudder angle called for is pro- 
portional to a combination of heading deviation and rate of change of 
heading, and which has a behaviour with a close resemblance to the 
automatic pilots used in practice. 
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Most modern steering engines are designed to move the rudder 
with an essentially constant speed, the rudder turning at that speed 
as long as a control signal is transmitted to the steering engine. A 
gyro pilot may be used to switch on this signal, calling for a correct- 
ing rudder at a certain small deviation from the desired heading, 
whereas a contact on the rudder may stop the motion at a suitable 
angle, the rudder remains in this position until the ship has swung 
over to the other side of its course and the rudder is then reversed. 
More often a follow up mechanism assures a "proportional control" 
of the rudder. Due to the finite rudder speed, these systems may be 
self exciting, however, and excessive oscillations may be built up. 

In order to overcome these difficulties the automatic pilot must 
be made to anticipate the motion of the ship, much in the way an 
experienced helmsman gives an auxiliary rudder. In practice this is 
accomplished by means of a feedback of rudder motion to the heading 
error detector, as in many commercial applications, sometimes also 
by adding some kind of rate of change of heading control, similar to 
the pitch rate component of submarine depth control systems. In 
effect, both these methods correspond to a "proportional plus rate 
control;" the first one often incorporates a. non-linear character to 
the system [ 16]. 

With these concepts in mind we are going to try to stabilize our 
unstable system, first by assuming that the autopilot is just a gain G, 
we represent the block diagram in figure 4 




FIGURE 4 



The open loop transfer function for this svstern is: 

\-\cs) _ <5 f< (_ 

Sit*- sr^Cf+sTt) Ct+sie) 

JKrom experience it has been determined that T ranges from 

ilr 

1 to 2, and a good choice is 1. 7 seconds, so the only parameter avail- 
able to make the system stable is G. 

Since there is one pole in the right hand plane, given by T^, there 
is no sense in using the Bode plot to determine stability values. 

1. Ro ot- Locus Analysis 

Getting the characteristic equation and arranging it in the 
proper form we obtain: 
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Since we have four poles there are four separate loci. 



There are three asymptotes due to 
n # of poles = 4 
q ir- of zeros = 1 
Center of asymptotes = 



n-q = 3 

Pj “ — - (3. 1 / / 
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Angle of asymptotes with respect to the real axis: 
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From a sketch made for this analysis we get the cross point 
at the jw axis = . 025. So G at the crossing point: 

^ ' |sl Uto.5" < jH5"<5.oo37| h+Odolf/ 
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2. Routh Criterion 

To find a better set of numbers we refer to the Routh 
criterion: ■ 
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For the stability limit in the s row, B must be equal to 



zero so; 



A (- 
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after manipulation: 



G - 85 G + 240 = 0 



from where 
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9 “ 3 at the crossing points 
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So we know that we have two crossing frequencies in the j\v 
2 

axis. In the s row 

A s Z + 1. 02 x 10" 5 G = 0 

replacing the values of G, and after all the mathematics we obtain: 
j .0225 

jw = 

J j • 15 

values that agree with those approximated from the sketch. 

In the figure 5, computer output of the subroutine Root Locus, 
we obtain the open loop characteristic of this system. As we can see 
the system can be considered as marginally stable due primarily to the 
long transient that can be predicted from the root locus graph. 

Using the DSL package program we simulate the system, 
first the system represented by the block diagram of figure 4, and 
after that the same system with the concept of ’’proportional plus 
rate control, 11 using rate of turn as feedback, as represented in 
figure 6 




FIGURE 6 
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In both cases, the disturbance was simulated as a step of 
amplitude 1, values of G were selected between 12 and 36, which are 
the ones with better time constant, copy of both computer programs 
is given at the end of the thesis. 

Figures 7-12 are computer outputs of the first condition, 
i. e. , system of figure 4, we can see as expected, the long time 
required to settle down due to the small time constant of the G values, 
but even so we can see in figure 9 for the value of G of 24. 2, com- 
paratively a short time for reaching steady state. 

Figures 13-18 are computer outp\ats of the second condition, 
i. e. , system of figure 6, we get better damped results, with a sig- 
nificant decrement in the time to reach steady state as exj^ected. 

Also in figure 15 we can see that comparatively this value of G gives 
us the best results. So we can state that the value of G = 24. 2 is a 
good value to stabilize the system and perhaps the best under this 
conditions. 
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IV. BASIC COURSE-KEEPING REQUIREMENTS 



It is well known that the course-keeping quality and the turning 
quality of a ship are contradictory to each other, and that it is im- 
possible to improve both of them without increasing rudder area [4]. 

The success of a rudder, designed for a particular ship, is 
measured by the degree to which the ship achieves a desired or anti- 
cipated course-keeping and course- changing ability. Since most ships 
spend most of their operating time moving ahead, it is assumed that 
maneuvering performance is specified only for the ahead direction. 

Course-Keeping ability may be specified quantitatively in terms 
of the stability index, of the characteristics of the Dieudonne Spiral 
Maneuver (see Appendix B), or of the range of rudder angles used 
to maintain a straight course. In reference [18], it is suggested 
that an attempt be made to design all ships for a stability index of 
zero or less, but it*is recognized that this may not be practicable for 
all ships . 

A major reason for the specification of a level of control-fixed 
stability for all ships is the fact that excessive controls -fixed 
instability leads to excessive use of the rudder to maintain a straight 

2 

This measure of Course-Keeping ability is frequently used in 
maneuvering in restricted waterways or in rough seas. 
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course. 



This in turn leads to increased wear on the rudder system, 



decreased ship speed, or increased fuel consumption. Controls- 
fixed instability also leads to increased difficulty in navigation in 
restricted waterways and in following seas and in avoiding collision 
with passing ships. 

A. UNDER STEADY CONDITIONS 
1 . Calm Water 

Among the factors important for the turning motion of ships, 
there are transverse forces and moments acting upon ship’s body 
itself. These may be approximately looked upon as motion in infinite 
fluid where image of ship is considered in case when V is not so large 
and the effect of waves by the ship itself with Froude Number below 
about 0.2 or 0.7 V L/2d [16] is negligible or when motion of ship 
is slow, and are the functions of V, ft , Q . This will be sufficient 
for turning motion of ordinary ship. 

The disturbances acting upon the system can be classified 
according to their influence on the behaviour of the system as follows: 

a. Disturbances that cause deviations from the set course. 

b. Disturbances which affect the steering characteristics of 

the ship. 

Wind and waves belong to the first group, to the second class 
of disturbances belong the loading of the ship, the depth of the water, 
etc. 
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In considering the calm water condition we arc assuming no 
disturbance at all, so, under these conditions we expect the ship to 
maintain its heading, without any action of the autopilot. 

2. Steady Wind 

Due to the wind pressure, the ship will be acted upon by a 
lift and drag which will cause the ship to drift at a drifting angle. By 
the drifting, the hull together with the rudder will create a certain 
hydrodynamic lift, drag and moment until the windward component of 
the hydrodynamic force balances with the wind drag. At this condition, 
the total summation of forces including ship { s thrust is zero. How- 
ever, the summation of the moment including moment due to the rudder 
at maximum helm angle will not necessarily be zero. 

If the summation of the moment is not zero, and is in a direc- 
tion to increase, the ship will have a tendency to turn windward. If 
the moment is in an inverse direction the ship will turn leeward. In 
both cases, the ship is supposed to be uncontrollable. If the summa- 
tion of the moment is zero within the range where rudder angle docs 
not exceed 30 degrees, the ship can keep its heading so that it is 
supposed to be controllable. 

B. ACCURACY NEED FOR COURSE- KEEPING 
1. Hoop Gain 

The most important parameter to be considered in considering 
automatic control is the Loop Gain, and its value depends directly on 
the transient response of the system. 
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For linear systems, the most commonly used correlation 

between frequency response and transient response is that correlation 

which exists between the height of the resonance peak M and the 

pw 

height of the peak overshoot of the step response M for a second- 
order system. A curve relating these is easily calculated [20j. 

Then, for a known second-order frequency response, the curve is 
entered with Mp W and the transient peak overshoot Mpj. is predicted 
exactly. When the system is known to be third-, fourth-, or higher- 
order, there is no such readily available correlation, but the second- 
order correlation may be used as an estimate. This gives a very 
accurate estimate if the higher-order response is dominated by one 
pair of complex roots, a less accurate estimate if dominance is not 
assured. 

We know that in order to make a system stable within some 
degree of accuracy we need to have complex roots, and a certain 
value of gain. In our study where we are operating with a giant of 
200, 000 Tons, , we can afford to have a certain degree of inaccuracy, 
i. e. , 1 to 2 degrees of heading error, without consequences, so we 
can go, if required, to have real roots governing our system. 

2. Loop Type Number 

Control systems, in general, are required to have certain 
operating characteristics which are determined quantitatively by the 
specific control problem, but which may be listed qualitatively the 
ones we are interested in as follows: 
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a. The system must be accurate in steady stale. 

b. The system must be stable. 

c. The system must regulate against disturbances. 

Each problem on feedback control has a different set of per- 
formance requirements that must be satisfied, but the requirement 
that is common to virtually all problems is a need for accuracy. 

What is meant by accuracy depends on the specific physical applica- 
tion, and no convenient broad definition is available, but the scope of 
accuracy requirements is readily illustrated. The class of control 
system we are interested in is the class of positioning system or 
servomechanisms. These systems are designed to change the output 
quantity as commanded by an input signal, and in addition are required 
to act as regulators in the presence of output disturbances. The 
primary consideration is usually that of keeping the error (difference 
between commanded output and actual output) less than a specified 
amount when the system is in steady state with no load disturbances. 

If static load disturbances are anticipated a separate specification 
for steady error is usually given, and for a suddenly applied disturb- 
ance a maximum instantaneous error man/ be specified. 

For positioning systems subject to command input signals, 
there are also specifications as to the permissible nature of the 
transient response, but these specifications usually refer to the 
permissible time duration of the transient and to the permissible 
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nature of any oscillations during the transient period. They seldom 
refer directly to the accuracy of the system, but the constraints 
which they place on the system design may make it very difficult to 
satisfy the accuracy specifications. 

Systems classified as positioning systems may be subjected 
to various types of input commands. These may be well defined and 
representable as simple mathematical functions or they may be rather 
complex functions of time. For well defined inputs the system usually 
(not always) reaches a finite steady state condition so that the accu- 
racy of the system is readily evaluated. For more complex inputs 
the system may not reach a clearly defined steady state condition, 
so that the accuracy must be evaluated in terms of bounds and a 
numerical value for M ste3.dy state error M may not exist. 

In this study emphasis is placed on system analysis and 
design for cases where the accuracy problem is defined in terms of 
deterministic input signals. 
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V. DESIGN OF COMPENSATION 



Having been able to obtain the best value for the gain of the 
autopilot to stabilize the unstable model, we see that the transient is 
still too large and has overshoots (due primarily to the hydrodynamic 
coefficients of the ship, which are given, as we may recall by the 
steering quality indices). 

With the concepts of chapter IV we are going to approach the 
design of the compensation using the Root- Locus method instead of 
the Bode Plot, since having a pole in the right hand plane the latter 
is hard to use. 

A. RATE GYRO COMPENSATOR 

Because of the shape of the Root Locus diagram we need to 
improve the transient response. By using a Rate Gyro we introduce 
an additional zero at 1/k, and this relocates the roots on the Root 
Locus plot at a different gain level. Pvate feedback is a very common 
means of increasing equivalent fiscous damping and of thus improving 
system transient response. Figures 19 and 20 are a block diagram 
representation of heading feedback, the original circuit and the 
equivalent respectively. 

The new characteristic equation is 

y & '1 C f . ±A±l = a 

SCt+ *7/ ) ( / 1 sTt.) ( / * sT ±) 
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FIGURE 19 




FIGURE 20 
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Figures 2l_26 are computer outputs of the subroutine Root-Locus. 
They give the different system behaviour for the values of k selected. 

As we can see from figures 24-26 in order to have a good time constant 
we need a large value of G, and from fig\xres 21_23 we can select 
very small values of G, which will give us real roots. 

B. LEAD TYPE FILTER COMPENSATOR 

Another type of compensation uses a filter with transfer function 
of the form: 

+ p 

having high pass characteristics . We are going to relocate, also, 
the roots on the Root-Locus diagram. This filter is just a cascaded 
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transfer function in the direct path between input and output, on our 



system. We are going to introduce one additional pole and zero. 

Figure 27 shows the block diagram with the filter in the system 




FIGURE 27 
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Figures 28-33 are computer outputs of the subroutine Root-Locus. 
They give the different system behaviour for the values of z and p 
selected. We can see from figures 31-33 that in order to have a 
good time constant we need a large value of G, while from figures 
28-30 we can select very small values of G, which will give us real 
roots . 

Having obtained the computer outputs for both types of compensa- 
tion, we can observe a similarity in the system behaviour for same 
roo t location. So we can predict the same transient for step inputs. 
However, we need to find an optimum value of k in the case of the 
Rate Gyro and of z for the Filter, in order to expand the little loop 
in figures 21-23 and 28-30. This optimum is the best value possible 
to obtain and it is a compromise between a. high gain and a good rudder 
operation, bearing in mind that for this to occur we must consider 
maximum deflection and rate for the rudder of this ship. Since for a 
quick response we need a high gain, but this could demand an exces- 
sive operation of the rudder and this could be bad for the following 
reasons : 

1. To get quick response we need faster rudder angle and a 
greater rudder angle. 

2. Faster responses may cause accelerations that are too 
sharp for personnel on ship. 
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So the next step will be to try to find the optimum value men- 
tioned above. Right off hand we can consider the system behaviours 
of figures 25, 26, 32 and 33, are very much similar to our original 
system without compensation, so we can disregard those four systems, 
and now we have four different systems for the Rate Gyro and Filter 
compensators respectively. 

We are going to simulate a disturbance that produces a rate of 
turn of .2 degrees /second, a value that has been found by experience 
can easily be reproduced on this tanker, requiring a correcting rudder 
angle of 10 to 20 degrees. 

With this external disturbance represented as a step, we examined 
our eight compensated systems. By trial and error method we check 

> 

for the optimum value required and we end up with a 1/k = . 041 for 
the Rate Gyro, and a z = . 041 for the Filter, which in essence are 
the same values we have for the case of 1/k = .04 and z = . 04 for the 
Rate Gyro and Filter respectively, so now remain these two systems 
to be studied for the best gain. In all the eight initial cases we ob- 
served a common characteristic, with minor variations, that for a 
small variation of heading, a sudden change in rudder occurs which 
causes a very quick response and we may suspect an excess in both 
rudder parameters. 

Figures 34-37 are computer outputs for the system compensated 
with the filter, we have chosen values of gain from the root-locus 




69 










w/ 1 • / l— 


RUDDLR 


RNGLf , THLTfl US T I ML 


K1 -9.42 


RGURVO 


PI GURU 34 






r,:0.20 



UN ITS/ INCH 



PLOT NO- L 



i 



70 



RUDDER RNGLE, THETfl US TIME Kl-lS-ci? 
RGURVQ FIGURE 35 



r ) r ) 

V 9 J 




VSCRL.r 0.40, LiNJTS/JHCH PLOT HQ. i 



71 



w ' 



RUDDLR ANGLE:, THLTA US TIME Kl-23 63 
RGURVO FIGURE 36 




72 



'/ 



RUDDER RNGLE, THEIR US TIME Kl-U-S 
RGURV0 FIGURE 37 




S'SCRLf -’O- 10 UNITS'' INCH PLOT NO- 1 



I 






73 



graph, for four different damping coefficients, 9.42, 15.87, 23.67 
and 3. 5 respectively. We can observe from the first three figures 
that a quick response is obtained and steady state is obtained in a 
very short time, but this could exceed the restrictions of the rudder. 

In figure 37 we see a moderate response and we do not reach the 
steady state at the end of the simulated time, which is only 5 minutes, 
as expected for a ship of characteristics such as this. We need to 
verify if there is any violation of the rudder restrictions. 

With the help of computer simulation, a program in DSL, which 
is given at the end of the thesis, is used to simulate limiters for both 
parameters. In order to introduce the limiting value of Rudder Rate, 
we need to replace, in the block diagram, the steering Gear transfer 
function for a circuit equivalent, given in figure 38: 




FIGURE 38 



where 



K 

3 



1/T 



E 



With the program above mentioned and the values of gain selected 
before we obtain the figures 39-42, computer outputs, where we see 
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in the first figures that even if the system is stabilized, there is a 



portion of the response where we exceed the rudder deflection, while 
in figure 42 we see a smooth settle down without any clipping. Will 
be a good idea to try values of gain within this small range of gain 
values . 

Figures 43-48 are computer outputs for the system above men- 
tioned with values of gain of 6. , 5. 3, 4. 6, 3. 5, 2. 5 and 2. We can 
observe in figure 44 for a gain value of 5. 3 we obtain a good settling 
time with the smaller steady state error of only . 8 of degree out of 
course which is excellent, but we have some rudder action even 
after settling time. In figure 45 for a gain value of 4. 6 we obtain a 
steady state error of 1. 03 degree out of course, with a settling time 
of about 2 more minutes than for the former case, but we get smaller 
rudder deflections. So we can conclude that this value of 4.6 will be 
-^alue for this system. 



the same disturbance, we obtain the computer outputs of figures 49-52, 
for values of gain selected from the Root-Locus graph, for four dif- 
ferent damping coefficient, 22.26, 67.28, 35.29 and 3.5 respectively, 
we observe almost the same results as in the case of the Filter, i. e. , 
the first three figures evidently violate the rudder restrictions and 
only in figure 51 for a gain value of 3. 5 gives us a good system 
operating under real conditions. Figures 53-58 are computer outputs 




Now considering the Rate Gyro having a zero at 1/k = . 04, with 
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for this system, for gain values of 6. 2, 5. 6, 5. , 4. 4, 3. 8, and 3. 2, 
following the same considerations as with the Filter. We can see 



from these figures that in all cases the system has increased its 
transient, being the best value, 6.2 as observed in figure 53. In the 



rudder action, so this value of gain is chosen as the best suited for 
this type of system. 

From comparison between these two systems we can select the 
one with Filter compensation, first, because of the best final results, 
i. e. , less amount of rudder required, better settling time, and 
secondly for the reliability of the system, since the Rate Gyro could 
not be suitable for marine use due to this overhaul interval. 



event of increasing the gain to 6.7 we start obtaining clipping in the 
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VI. CONCLUSIONS 



A. RESULTS 

An optimal value of gain can always be obtained using the Root- 
Locus method for any given plant. Both negative and positive feed- 
back can be used to achieve the desired results but only negative 
position feedback should be used to obtain a stable system. 

With the concepts of Chapter IV, the requirements for a basic 
Cours e-Keex>ing, and our best value for the system selected as the 
best, we verify with a low order disturbance represented by an initial 
condition of 0 , that the actual heading in the steady state is equal 
to the reference heading or which is the same the steady state error 
is equal to zero. We already have observed that for a load disturbance 
rex^resented by the step input we obtain a small steady state error. 

So we have reasons to believe that we are working with a type zero 
system and, if it is so, do we need to change the loop type number? 

By doing an analysis of the feedback loop we can easily verify 
that indeed we have a type zero system, and using the final value 
theorem for this linear system, the evaluation of steady state error 
for the deterministic input is easily accomplished. In this case 
where steady state error is finite its magnitude is determined by the 
reciprocal of the gain transfer function for the direct path between 
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input, and output, this magnitude is called an error constant, also 
called the DC gain or zero frequency gain, because it is the numerical 
value obtained for the transfer function by deleting the s factors and 
substituting s = 0. Following this reasoning and with the final results 
obtained we see that this meets our requirements for the stable 
operation of the system with the degree of accuracy that can be con- 
sidered as satisfactory. 

The purpose of this thesis, which was to study the basis of the 
Course-Keeping with Automatic Control, is reached at this stage. 

It is good to point out that further studies could be made by studying 

V 

the sea effects, how automatic steering and the yawing of ships is 
affected in rough seas. 
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LIST OF SYMBOLS 



draught of a ship, mean of fore and aft if no remark. 



moment of inertia of a ship about a vertical axis through 
her centre of gravity. 




moment of inertia of additional mass about a vertical axis 
through the centre of lateral additional mass. 



J; 



zz 






length of a ship, between perpendiculars, if no remark. 



mass of a ship. 

w/t lU 

longitudinal additional mass. 



lateral additional mass. 



w //t L ‘ d 

hydrodynamic moment about a vertical axis through the C. G. 
of a ship, positive for standard turning moment. 

//// L Z d ^ 1 

d"7is 



yaw angular velocity, positive to starboard, identical with 0 
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r/(V/L) or L/R 



R turning radius, steady or instantaneous 

V ship speed in meters/second 

Y lateral component of hydrodynamic force, positive to star- 
board. 

Y ' 

Yjj zY'kn 

Y 1 DV /c> -r' 

Y' s bV'/ts 

angle of drift, positive to port. 

£ angle of helm, positive to starboard. 

Q turning angular velocity, positive to starboard, identical to r. 
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APPENDIX B 



THE DIEUDONNE SPIRAL MANEUVER 

The spiral maneuver consists of the following: 

1. The ship is ''steadied 1 ' on a straight course at a preselected 
speed and held on this course and speed for about 1 minute. Once a 
steady speed is established, the power plant controls are not manip- 
ulated for the duration of the maneuver. 

2. After about 1 minute, the rudder is turned to an angle, £ , 

of about 15 degrees and held until the rate of change of yaw angle 
maintains a constant value for about 1 minute. 

3. The rudder angle is then decreased by a small amount 
(about 5 degrees) and held fixed again until a new value of <3 is 
achieved and is constant for several minutes. 

4. The foregoing procedure is repeated for different rudder 
angles changed by small increments from, say, large starboard 
values to large port values and back again to large starboard values. 

The numerical measures obtained from the preceding spiral 
maneuver are the steady yawing rates as a function of rudder angle. 

A plot of these values is indicative of the stability characteristics 
of a ship. For example, if the plot is a single line going from star- 
board rudder to port and back again, as shown for ship A in figure 1-B, 
the ship possesses controls-fixed, straight-line stability; that is, it 
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STAC , tft SHI? A 
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has a negative stability index. If, however, the plot consists of two 
branches joined together to form a "hysteresis 11 loop, as shown in 
ship B of figure 1-B, the ship is unstable; that is, it has a positive 
stability index. In addition, the height and width of the loop are 
numerical measures of the degree of instability; the larger the loop 
the more unstable the ship. The slope of the yaw- rate curve at zero 
rudder angle is a measure of the degree of stability or instability. 
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C THIS F R CG PAP WILL FLCT Tt-E SYSTEM RESPONSE TO A STEP 

C CISTUREANCE* THE ELCCK DIAGRAM OP THE SYSTEM IS THAT 

C REPRESENTED IN FIGURE*. THE VALUES OF GAIN SELECTED 

C' A RE THE ONES WITH BETTER TIME CONSTANT FROM THE ROOT 

C LCCLS GRAPH OF FIGURE 5 

PA RAM IM=C. , IN2 = 1. ,K1=12. ,K2=-.0434 

CCTPL F INT IM = £CC . ,DELT=.8 ,DELS = . 8 

INTCEF N CM, N PLOT 

CONST NP LC T = 1 

DERIVATIVE 

EPPCF= INI-THETA 
CELTAR=Kl^ERROR 
CELT A=REALPL( 0. , l.?,DELTAR) 

ANPL I=K2*CELTA 

CCMF=LEDLAG(0.,20.,-269.3,AMFLI) 

CCR R E ( - R E ALPL ( 0 . , 9 . 3 , COMP ) 

THETA 1= IN2+CCRREC 
TFETA = INTGRL(0. , Th ETA1 ) 

S A V F L E 

CALL CFNG(1,1, TIME, DELTA) 

CALL DRUGU, 2, TIME, THETA) 

TERMINAL 

CALL ENCPWNPLCT) 

END 
ST C F 
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C THIS PROGRAM WILL PLCT THE SYSTEM RESPONSE TO A STEP 

C DISTURBANCE. THE BLOCK DIAGRAM OF THE SYSTEM IS THAT 

C REPRESENTED IN FIGURE 6. THE VALUES OF GAIN SELECTED 

C ARE THE ONES KITH, BETTER TIME CONSTANT FROM THE ROOT 

C LCCLS GRAPH OF FIGURE 5 

PAR AM I M = 0 . » I N 2 = 1 . »K1 = 12. ,K2=-.C434 

CCTFL PINT IM = 800.fDELT=.8»DELS=.8 

INTGEP NLMtNFLCT 

CONST N F L C T = 1 

DER IV AT 1 V E 

EFFCR = I NI--THETA-THETA1 

DELTA R = K 1* ERROR 

CELT A=REALPL(0., 1.7, DELTAR) 

AMPL I =K2«CELTA 

CCKP=LEDLAG( C. ?2C. »-269.3.AMPLI ) 
CCRFEC=REALPL(0.»g.3,CCMP) 

TFETA1=1N2+CCRREC 
T F ETA= INTGRL ( 0 . t Th ETA 1 ) 

S A M F L E 

CALL DRWG ( I f 1 ,TI ME f DELTA) 

CALL CRWG( 1 * 2, TIME, THETA) 

T E F M I f\ A L 

CALL ENDRMNPLQT) 

ENC 

STOP 
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C THIS PROGRAM HILL PLOT THE SYSTEM RESPONSE TO A STEP 

C DISTURBANCE. THE ELOCK DIAGRAM OF THE SYSTEM IS THAT 

C REPRESENTED IN FIGURE 34. THE VALUES OF GAIN SELECTED 

C ARE TFE ONES WITH THE DAMPING COEFFICIENTS SELECTED 

C FRCM TFE RCOT-I.OCUS CIAGRAMS OF FIGURES 23, 25, 30, 31 

C AND UHH THE LIMITERS FOR TFE PUCDER PARAMETERS USING 

C TFE ECU I VALENT CIRCUIT OF FIGURE 43. 

// EXEC DSL 
//DSL. INPET CD * 

PA RAM IM = 0. ,IN2=«0033,K1=9.42,K2=— .C434,K3=.588, . . . 

F L=~ .C4 04 , P2=.0404fP3=-.52,P4= . 52, K4=0. 

CCNTFL FINTiM=80C.,DELT=.8,DELS=.8 
CCNS7 $hl=G, SW2=0.SW3=C» FLAG=O f NPLOT=4 
INTGER FLAG , Sh 1 , SW2#SW2 » NPLQT 
DERIVATIVE 

EFFO«5*INI-THETA-ThETAI*K4 

FILT5R=LEDLAG(0.»25.,2.5,ERRQR) 

DELREF=FILTEft-Kl 
ALFF A=CELREF-DELTAC 

PRCCED ALFA»RUCCtR( ALPHA , FLAG, SW1 t SW2,SW3, DELT AC » P4 ) 

ALFA =A L PhA 

IF ( A 8 S { CELTAC) .G E . P4 ) SH 1= I 
IFiSNl.EC.l) ALFA = 0 . 

IFCSH1 .EC.O) S W 3 = 0 
IF { ALFFA J 10,6,6 
10 SR 2= I 
GC TC 1 



6 $ H 2 = C 

7 IF ( ( ( FIAG.NE.SH2) .OR. (SW3.EQ.1 H .AND. ( SW1 . EQ . DJGO T08 
GO TO 9 

fi ALFA=ALPFA 
SHE = 1 
5 SH 1= C 
F L A G = S H 2 

ENDFFC 

CELDCR= K3-ALFA 
C ELDOT = L I P IT ( PL, P2, DEL COR) 

CELTA1=IKTGRLC0. , CEL DOT l 
CELT AC *LIPIT(P3,F4, DELT Al) 

AMFL I = K2 !> CELT AC 

CCMP=LECLAG(0. ,20. ,-269 .3, AM PL I) 

CCRREC=REALPL( 0. ,9 .3 ,COMP» 

TF ET A l=IN2 cC0RREC 
TFETA=INTGRL(0,TFETAiI 

SAMPLE 

CALL CRHG( 1, I, TIME , D E L T AC ) 

CALL C RWG ( l ,-2, T IM E , TFETA ) 

TERMINAL 

CALL ENDRMNPLOT) 

GC TC (I, 2, 3), N PLOT 
2 K 1 = 1 5 . £ 7 
CALL RERUN 
RETURN 
2 Hi = 23. 6 7 
CALL RERUN! 

RETURN 



1 K 1 = 3 . 5 
CALL RERUN 

END 

STOP 

//PLCT.SYSIN DD * 

RUDDER ANGLE, TFETA VS TIME 
AGLAYC FIGURE 39 

RUDDER ANGLE, TFETA VS TIME 
AGLAYC FIGURE 40 

RUCDER ANGLE, TFETA VS TIME 
AGLAYC FIGURE 41 

RUDCER ANGLE, TFETA VS TIME 
AGLAYC FIGURE 42 
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